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ON  THE  COMPUTATIONAL  PHOCEOTRES  FOR  TIRING  AND  BOMBING  TABLES 


ABSTRACT 

This  report  describes  the  methods  and  means  of  computing  ground-to- 
ground  firing  tables  and  air-to-ground  bombing  tables  at  the  Computing 
Laboratory,  BRL.  The  description  includes  all  steps  in  the  transformation 
from  measured  data  into  printed  tables.  It  includes  a detailed  discussion 
of  the  effect  of  the  use  of  high  speed  digital  computers  in  various  portions 
of  the  computations,  and  discusses  further  the  expected  effect  of  further 
"mechanization". 
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Intentionally  deft  blank. 


INTRODUCTION 


This  report  is  primarily  concerned  with  the  description  of  the  computa- 
tional procedures  used  at  present  at  the  Aberdeen  Proving  Ground  to  produce  firing 
and  bombing  tables.  The  data  gathering  methods  and  the  field  requirements  on  the 
tables  will  be  mentioned  only  insofar  as  they  affect  these  procedures. 

These  present  procedures  were  developed  continuously  over  a number  of  decades. 

At  the  beginning  of  the  first  world  war,  except  for  the  case  of  mortar  fire,  the' 
only  trajectories  being  computed  in  this  country  were  for  direct,  low  angle  fire, 
and  the  Siacci  methods  were  used  (see  Ingalls”  Tables  - Artillery  Circular  M). 

With  the  introduction  of  anti-aircraft  guns,  it  became  apparent  that  this  was 
insufficient.  At  Sandy  Hook  and  at  Aberdeen,  Lt.  Philip  Alger  began  the 
computation  of  AA  trajectories  by  the  French  Short  Arc  methods  in  1917-  Early 
in  1918,  Capt.  F.  W»  Loomis,  who  had  visited  the  French  and  British  ballisticians, 
introduced  at  Aberdeen,  Fowler's  method  of  numerical  integration.  Prom  then  on 
the  contributions  of  well  known  mathematicians  to  the  development  of  the  present 
conputational  methods  were  continual,  beginning  with  Major  Veblen’s  Range  Firing 
and  Computing  Section  at  Aberdeen,  and  followed  shortly  thereafter  by  Major  Moulton’s 
group  at  the  Office  of  Chief  of  Ordnance  in  Washington.  Veblen’s  group  included 
Bliss,  Wiener,  Gronwall,  Franklin,  Gill,  Widder,  and  H.  H.  Mitchell.  Moulton’s 
group  included  Ritt,  A.  A.  Bennett,  Milne,  D„  Jackson,  Ford,  Alexander,  Buck, 

Walsh,  Vandiver,  Dresden,  and  Dantzig.  In  the  past  decade,  further  development 
of  these  procedures,  including  the  development  of  bombing  table  computations, 
paralleled  the  growth  and  emergence  of  the  Computing  Laboratory  as  a separate 
unit  under  A.  A.  Bennett,  T.  E.  Sterne,  L.  S.  Dederick,  and  W„  W.  Leutert. 

However,  since  the  developments  and  refinements  in  the  overall  procedures  appeared 
at  different  times,  until  now  no  complete  picture  has  been  given.  To  give  such 
a picture  of  the  procedures  used  to  obtain  two  fairly  typical  kinds  of  tables  is 
one  of  the  purposes  of  this  report.' 

This,  of  course,  does  not  mean'  that  developments  and  refinements  in  these 
procedures  will  not  continue.  On  the  contrary,  t he  pressure  of  the  development 
of  new  weapons,  the  discovery  and  improvement  of  measuring  instruments  and  tech- 
niques, the  advancement  of  the  physical  theories,  the  evolution  of  new  mathemati- 
cal methods,  and  the  construction  and  procurement  of  more  powerful  computing 
devices  will  force  these  computational  procedures  to  change.  . Indeed,  another  pur~ 
pose  of  this  report  is  to  serve  as  a reference  for  those  who  are  in  a position 
to  encourage  and  implement  such  future  improvements  of  the  procedures  as  result 
from  the  above. 

In  addition,  this  report  serves  as  a more  detailed  description  of  the 
execution  of  an  essential  mission  of  the  Computing  Laboratory  and  acquaints  new 
personnel  with  the  details  of  this  mission,. 

From  the  statement  of  the  above  three  purposes  for  this  report  it  is  clear 
that  it  is  not  intended  to  establish  or  advocate  a standard  operating  procedure 
for  the  production  of  Ballistic  tables}  nor  is  it  intended  to  recommend  per- 
petuation of  present  practices. 

By  the  "computational  procedures"  to  have  been  talking  about  we  mean  the 
complete  programs  for  transforming  the  measured  data  into  tables;  these  programs 
involve  the  use  of  many  sub-procedures  including  hand  computations  and  various 
computations  by  possibly  -the  same  or  different  high  speed  digital  computing 
machines.  Very  many  of  the  overall  procedures  have  been  recently  mechanized, 
and  more  mecnanization  is  certainly  possible.  A further  purpose  of  this 
report  is  to  point  out  where  Ihrther  mechanization  is  possible.*  However,  let 


us  remark  here  that  the  chain  of  a completely  mechanized  program  should 
be  broken  at  certain  critical  points*  if  only  to  examine  for  the  effects 
of  bad  data  which  may  need  rejecting,  "With  machines  having  available 
much  larger  memories  than  the  present  ones  nevertheless  one  can  see  that 
even  this  examination  for  bad  data  could  be  mechanised. 

The  computational  procedures  we  will  describe  are  those  of  the 
classical  problems  of  ground  to  ground  firing  tables  where  the  shell '3 
trajectories  are  not  tracked  on  the  range  and  tables  for  bombing  of 
fixed  targets  from  aircraft  where  the  bomb's  trajectories  are  tracked 0 
There  are,  of  course,  other  ballistic  tables  j but  as  will  be  pointed 
out  in  the  concluding  remarks*  many  of  the  features  in  their  production 
appear  in  the  production  of  the  above  two  types  , The  terms  firing 
and  bombing  tables  without  further  qualification  will  in  this  report 
always  refer  to  these  two  types. 

Let  it  suffice  for  the  moment  to  say  that  when  the  block  flow 
charts  of  the  computational  procedures  for  firing  and  bombing  tables 
(see  figures  1 and  2)  are  compareds  as  is  to  be  expecteda  striking 
similarities  in  structure  are  apparent , There  are  also  striking 
differences*  one  of  the  most  readily  observable  ones  being  the  ability 
to  adjust  the  drag  coefficient  for  each  bombing  table  because  the 
missile  is  tracked*  while  for  firing  tables  one  can  only  assume  a 
drag  coefficient  and  adjust  the  ballistic  coefficient  because  at  pre- 
sent the  missile  is  not  tracked,  This  explains  why  in  the  reductions 
at  present  the  bombing  tables  branch  requires  more  personnel  than  the 
firing  tables  branch.  The  use  of  accelerometers  in  range-dropped  bombs 
to  determine  the  drag  on  the  bomb  and  the  possible  use  of  Doppler 
techniques  to  track  gunfire  may  make  the  work  of  reduction  per  missile 
about  the  same  for  both  branches.  However*  it  should  be  remarked  that, 
since  the  present  firing  tables  are  more  than  accurate  enough  in  view 
of  varying  meteorological  conditions  and  the  lack  of  sharp  control  over 
the  muzzle  velocities*  tracking  gunfire  missiles  is  an  unnecessary  expense. 

Acknowledgement  of  the  authors5  gratitude  are  due  to  J,  Prevas*  H* 
Reed  and  others  of  the  Firing  Tables  Branch,  and  to  E,  H„  Martin,  Mo 
Field  and  others  of  the  Bombing  Tables  Branch  for  Extensive  time  and  per- 
tinent details.  The  authors  also  thank  L»  Butler,  E„  Gersten  and  others 
of  the  Coding  Sections  for  details  of  the  programs  and  their  running 
times o They  are  also  grateful  to  A,  A,  Bennett  and  I,  E,  Sterne  for 
their  critical  reviews  and  advice, 

II  - THE  TRAJECTORY  EQUATIONS 

The  physical  assumptions  behind  the  present  computational  procedures 
producing  firing 'and  bombing  tables  are  that  the  shell  or  bomb  are  essen-= 
tially  particles.  For  the  tables  discussed  in  this  report  the  only  con- 
cession in  the  equations  used  in  the  computation  to  the  fact  that  the 
projectile  is  not  really  a particle  is  that  there  is  a drag  force.  This 
force  per  unit  mass  is  given  by  a resistance  function  one  of  whose  factors 
is  a drag  coefficient.  Except  in  air  to  air  trajectories  * no  com- 
putationally feasible  theory  for  spin-and  finest ablized  shell  or  bombs 
with  a general  vaw  is  in  use  now.  Small  yaw  theories  have  existed  for 
some  time  [3] * But  because  of  their  computational  difficulties 

the  effects  of  yaw  on  the  trajectory 


Such  numbers  refer  to  the  bibliography 


have  only  been  computed  in  some  special  cases o (We  remark  that  at 
present  studies  to  incorporate  the  effects  of  yaw  in  more  general 
cases  are  continuing  flQ ) « 

The  equations  for  a particle  trajectory  with  drag  and  including 
(see  remarks  on  omissions  below)  the  effects  of  wind  and  the  Coriolis 
force  due  to  the  earth ss  rotation  are 

V **  = E (x  • w^)  + X^  y, 

# . <i  9 _ © O 

(1)  y 

OO  _ / O \ « 8 

z a - E (z  “ vz>  + \^y  *■  XgX 

where  x is  distance  down  range,  y is  vertical  distance,  z is  horizontal 
distance  to  the  right  (looking  down  range),  (\^9  - Xg,  - \^)  are  twice 

the  values  of  the  components  of  the  earth ffs  angular  velocity  in  the 
(x,  y,  z)  system,  E is  the  resistance  function,  and  g is  the  accele- 
ration due  to  gravity.  The  vertical  wind  component,  and  the  Coriolis 
contributions  involving  z are  omitted,  being  very  small  compared  to  the 
remaining  terms  in  the  equations  (l).  (In  the  bombing  reductions  the 
range  grid  coordinates,  which  differ  only  by  a slight  rotation  from 
the  above,  are  used.) 

For,  firing  tables  g is  treated  as  a constant  except  for  trajectories 
with  very  high  summits  (of  the  order  of  five  miles).  For  bombing  tables 

(2)  g gQ  (l  - “?) 

where  gQ  is  a constant  and  r is  the  earth fs  radius.  The  resistance 

function  E is  given  by 

, #>(y)  u Ktj(M) 

(3)  E - — — 

C 

where  p(y)  is  the  air  density  at  the  height  y,  u is  the  projectile *s 
speed  relative  to  the  air,Kjj(M)  Is  the  drag  coefficient,  an  empirical 

function  of  the  Mach  number,  M ( ®u/a),  of  the  projectile,  and  C is 
the  ballistic  coefficient  which  is  given  by 

(it)  C = 

id 

where  m is  the  mass  of  the’ projectile,  d3 is  its'fcaliber  and  i is  the 
form  factor  (a  constant  near  unity). 
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For  bombing  tables  it  is  customary  to  write 


($)  E = y p(y)  u Kp 

where  y is  the  reciprocal  of  the  ballistic  coefficient*  (Sometimes  p(y) 
is  replaced  by  H(y),  the  ratio  of  the  air  density  at  the  height  y to 
that  at  sea  level  and  then  y is  redefined  as  the  quotient  of  the  density 
at  sea  level  by  the  ballistic  coefficient*  This  is  done  to  replace 
handling  the  large  number  of  the  very  small  values  obtained  for  the  den- 
sity, p(y)s  by  handling  the  same  number  of  moderate  values  obtained  for 
the  relative  dens it  y3  H(y)0) 

For  firing  tables  two  forms  of  E are  used.  In  the  tradition  of 
the  work  of  the  Gavre  Commission  in  the  late  nineteenth  century  and  with 
later  modifications  to  include  the  dependence  of  the  drag  on  the  speed 
of  sound  E has  been  written  as 

(a/a  ) pfr)  G(u  /a2) 

(6)  E . 2 — — — 

where  a/a  is  the  ratio  of  the  speed  of  sound  a at  the  height  y to 
s 

that  a at  sea  level  and  G is  called  the  drag  function,  an  empirical 
s 

function  originally  considered  to  depend  only  on  u and  now  tabulated 
as  a function  of  (u/a  )«  When  used  in  computation  u2  is  multiplied 
2 ^ 

by  (a Jb.)  and  this  product  used  as  the  independent  variable  for  entry 

into  the  G - table,  the  form  more  in  line  with  the  current  aerodynamic 
usage  is  that  of  (3)* 

The  trajectory  equations  (1)  are  solved  numerically  for  x,  j9  and 
z by  the  modified  Euler  method  of  numerical  integration  with  one 
iteration  (sometimes  called  the  Heun  method) o The  general  procedure 
is  as  follows*  First  the  system  of  equation  to  be  integrated  is 
written  as  a system  of  first  order  equations  of  the  type 

(7)  V(t)  = F (V,  t) 

where  V and  F are  vectors*  Then  a first  approximation  to  V(t+  ^kt)  is 
given  by 

(0)  T (t*<4  t)  = V(t)  + F (V,  tM  t* 

the  next  and  final  approximation  is  given  by 

(9)  V(t*<4t)  = V(t)  + [F(V,t)  + F(V,  t*  t)]^t/2* 

In  the  case  of  the  so-called  reduction  trajectories,  described  in  sec- 
tion IH,  V and  F have  six  components,  namely 


V = x 

1 


V7 


fl"V2 


V2  f2  - - E [V2  - ». 


x 3 


V\ 


"Xi  ^ 


V,  - 7 


- 6 “ X1V2: 


7$  fc  z f g * V6 

V6  » i f6  - E [t6  - ..Or,,]  *yu  *x2t2 


In  the  case  of  the  so-called  normal  trajectories,  V and  F have  four 
components,  the  last  two  components  being  deleted  as  well  as  the  terma 
in  the  f*s  which  involve  the  w5s  and  the  \rs.  In  each  case  A t is 
varied  with  the  speed,,  This  is  necessary  because  the  usually  en- 
countered Kp  vs  o M curve  is  slowly  varying  up  to  about  M ■ ,?  or  so| 

then  it  rapidly  takes  on  a relatively  large  positive  curvature  and 
then  rises  sharply  somewhere  between  M = ,96  or  so  and  M = 1,2  or  so; 
then  the  curve  takes  on  a large  negative  curvature  and  finally  settles 
down  as  a slowly  varying  function  again.  The  range  of  the  argument 
of  Kp  is  divided  into  five  speed  intervals  characterized  by  the  above  • 

behavior.  At  each  step  in  the  integration  a discrimination  is  per- 
formed to  determine  in  which  of  these  five  intervals  the  speed  of  the 
projectile  lies.  For  all  the  steps  for  which  the  projectile  is  in  a 
particular  speed  interval  ^t  is  a constant.  These  five  constant 
values  of  A t have  been  determined  by  experience  with  many  K^'s.® 

We  add  here  that  very  extensive  discussions  of  the  theory  and  of 
the  hand  computation  of  firing  tables  are  given  in  the  books  by  Moulton 
m9  Bliss  [Vj  Hayes  [ll  , Kelley , McShane  and  Reno  (5)  and  the  roeroe- 
ographed  notes  oy  Dederick  CO?  and  the  report  by  Hitchcock  |_13j  . 


HI  - FIRING  TABLES 

A,  General  Remarks 

Firing  tables,  generally  speaking,  are  tables  of  data  for  use  in 
aiming  and  timing  bullets,  shell,  rockets,  or  other  projectiles®  We 
will  limit  most  of  our  discussion  of  firing  tables  to  tables  used  in 
ground  gun  fire.  In  addition  to  range  versus  elevation,  drift,  and 
time  of  flight,  among  the  most  useful  quantities  in  such  tables  are 
probable  errors  in  range  and  deflection,  range  effects  of  increase  in 
muzzle  velocity  and  in  rear  wind,  change  in  range  for  1 mil  change  in 
elevation,  and  the  like.  These  data  are  needed  by  artillerymen,  manu- 
facturers of  munitions  and  of  aiming  equipment,  and  others*  These 
tables  must  be  produced  by 
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ac  executing  a series  of  firings  and  making  relevant  measurements 
including  measurements  of  the  various  non-standard  conditions  at 
firing  timep 

b0  cbmputing  a series  of  pertinent  trajectories  to  estimate  the 
ballistic  coefficient  C and  the  effects  of  variation  in  C (These 
are  the  so-called  "reduction  trajectories'*) a 

Co  computing  a series  of  trajectories  under  "normal"  conditions^ 
i*e,5  with  v = w = X,  ®X«=X.j“0  and  a standard  atmosphere 

for  H(y)  (these  are  the  so-called  "normal  trajectories" ),  and  then 
the  effects  of  variation  in  weighty  muzzle  velocity,  atmospheric 
conditions,  and  angle  of  elevation  on  the  range s and 

d0  computing  a correction  formula  for  time  of  flight  and  a for- 
mula for  drift,  and  finding  probable  errors  in  range, 

A block  flow  chart  describing  these  steps  in  further  detail  may- 
be found  in  Fig,  1, 

Let  us  now  consider  what  is  done  in  the  main  blocks  of  this  flow 
chart,  and  how  the  results  are  obtained, 

B,  Collection  of  Range  Data? 

At  the  firing  range  the  position  of  the  trunnion  of  the  gun  is 
measured  by  a transit.  The  azimuth  and  elevation  of  the  gun  are  set 
and  measured  using  a transit  for  the  azimuth  and  a muzzle  clinometer 
for  elevation.  The  projectile  is  then  weighed,  magnetized,  and  placed 
in  the  gun.  The  time  of  firing  is  recorded.  The  initial  speed  is 
determined  by  the  use  of  a chronograph  and  two  parallel  coils  suspended 
ahead  of  the  muzzle  of  the  gun  in  such  a manner  that  the  axes  of  the 
coils  and  the  muzzle  of  the  gun  are  in  line,  (This  last  set  up  is 
called  a "cage").  The  magnetized  projectile  induces  currents  in  each 
of  these  coils  at  slightly  different  times  which  are  measured  by  the 
chronographj  this  time  difference  together  with  the  distance  between 
the  coils  is  sufficient  to  determine  the  average  speed  through  the 
cage.  Since  because  of  blast  the  true  speed  cannot  be  measured 
exactly  at  the  muzzle,  an  effective  muzzle  speed  is  obtained  by  an 
extrapolation  from  the  average  speed  through  the  cage.  The  impact 
point  which  is  usually  on  water,  is  observed  from  three  or  four  towers 
using  a theodolite,  and  the  time  of  flight  is  determined  by  stop 
watch,  (Whenever  the  projectile  is  fuzed  to  burst  in  the  air,  the 
time  of  flight  to  the  burst  point  is  determined  by  a chronograph,) 

These  data  are  gathered  for  various  muzzle  speeds  and  initial  angles 
of  elevation,  there  usually  being  approximately  ten  rounds  per  group 
with  five  or  six  different  elevations  for  each  muzzle  speed  and  fron 
one  to  eight  different  muzzle  speeds, 

While  the  firings  are  going  on  meteorological  data  (often  briefly- 
called  metro  data)  are  gathered  hourly  from  radiosondes  in  ascending 
balloons  which  are  tracked  by  direction  finders  and  theodolites. 


These  data  include  range  wind  w , cross  wind  w s air  temperature 

X z 

T,  and  humidity,  all  as  functions  of  altitude  yc  The  height  of  tides, 
which  is  necessary  to  get  the  height  of  impact,  (most  of  firings  are 
over  water)  is  obtained  at  the  same  time  from  a gauge  aticko  (Since 
the  purpose  of  this  report  is  to  present  the  logical  aspects  of  the 
computational  procedure,  the  reader  interested  in  further  details  of 
the  data-gathering  techniques  if  referred  to  the  many  pertinent  BRL 
reports  and  such  comprehensive  books  as  that  of  Crana,|lll  «.) 

These  raw  data  must  then  be  processed  by  preliminary  reductions 
to  determine  the  parameters,,  functional  as  well  as  numerical,  which 
enter  into  the  reduction  trajectories,  for  the  computation  of  the 
ballistic  coefficient,  C» 

C , Preliminary  Reductions  2 

The  numerical  parameters  are  the  mass  of  the.  projectile  and  initial 
conditions  required  for  solving  the  trajectories  (l),  namely,  muzzle 
coordinates,  x , yQx  zq,  muzzle  speed,  V , elevation,  0,  and  muzzle 

velocity  components,  is  s y , which  except  for  0 are  obtained  by  simple 

calculations  on  the  raw  data,  however,  the  0 obtained  from  the  muzzle 
clinometer  measurements  is  corrected  by  comparing  ranges  of  a few  com- 
puted low  angle  trajectories  with  observed  ranges * (This  correction 
is  called  jtmqps  it  covers  the  result  of  many  effects  of  the  vertical 
motion  of  the  gun  during  firing* ) 

The  functional  parameters  required  are  the  drag  coefficient  Kjj(M), 

or  the  drag  function  G,  the  air  density  p(y),  the  speed  of  sound  a(y), 
the  absolute  temperature  T(y),  and  the  wind  components  w (y)  and  v (y)* 

The  temperature  T(y)  is  used  to  obtain  a(y)  through  the  relation 


where  the  subscript  s denotes  that  the  quantity  modified  by  the  sub- 
script is  evaluated  at  sea  level  under  standard  conditions*  It  is  also 
used  with  the  humidity  and  pressure  to  obtain  the  density  p(y)» 

Since  the  reduction  trajectories  are  to  be  trajectories  represent- 
ative of  a group  of  firings  under  similar  conditions,  those  numerical 
and  functional  parameters  which  vary  from  shot  to  shot  are  averaged. 

The  various  averagings  are  performed  as  follows* 

The  mean  muzzle  velocity,  v , for  a group  is  computed  as  the  arith- 
metic mean  of  the  observed  values  for  all  the  firings  of  the  group. 

Since  the  elevations  are  reset  after  each  firing,  and  the  jump  adjust- 
ment is  uniform,  this  implies  that  the  muzzle  velocitity  components  for 

V " are  also  the  arithmetic  means  of  the  components  of  the  observed 
0 

muzzle  velocities,  corrected  by  jump* 


The  impact  points  for  the  firings  at  the  Aberdeen  Proving  Ground 
are  usually  on  water , Tne  height  of  tides  for  the  group  is  taken  as 
the  height  at  the  mean  time  of  fire,  This  is  used  to  determine  the 
mean  height  of  impact  for  the  group. 

The  terminal  values  of  the  time  of  flight,  t^*  the  range*  x^,  and 

the  deflection*  as  well  as  the  mass  of  the  projectile  for  our  aver- 

age  group  trajectory  are  also  determined  by  taking  arithmetic  means. 

The  group  time  of  flight*  tffl>  are  used  to  find  corrections  for  the  times 

of  flight  to  be  entered  in  the  tables.  The  group  ranges*  x^*  are. used 

in  the  criteria  establishing  the  form  factor*  i*  during  the  running  of 
the  reduction  trajectories.  They*  as  well  as  the  group  deflections* 
a * and  the  sampling  statistics  to  be  described  shortly  are  used  in 

the  computation  of  probable  errors  in  range  and  deflection. 

The  averaging  procedure  for  the  metro  data  (namely  T*  v * w * and  p) 

X 2j 

is  slightly  more  complicated.  First  the  mean  times  of  firing  for  those 
portions  of  the  group  of  rounds  between  successive  metro  ascensions  is 
foundj  then  the  metro  data  are  interpolated  for  equal  intervals  in  alti- 
tude at  each  ascensions  then  a further  interpolation  is  made  to  the 
mean  times  for  these  equal  intervals  in  altitudes  and  finally  the  metro 
data  at  equal  y for  the  various  mean  times  are  weighted  by -the  number 
of  rounds  between  the  various  successive  ascensions  to  obtain  the 
mean  metro  data  for  the  whole  group. 

With  the  completion  of  hll  these  reductions  of  the  raw  data*  which 
are*  incidentally*  being  done  by  hand*  a drag  coefficient  or  drag  func- 
tion* Kp  or  G is  chosen  from  among  many  available.  This  choice  is 

based  upon  experience  in  the  past  with  projectiles  of  all  shapes  and 
sizes. 


The  final  reduction  made  before  a first  estimate  of  the  ballistic 
coefficient  C is  made  and  the  reduction  trajectories  run  is  the  analysis 
of  the  observed  probable  errors  of  range  and  deflection  made  to  minimize 
possible  trends  in  wind  and  velocity.  This  is  done  as  follows.  First 
the  usual  sanpling  estimate  of  the  variance  is  obtained  s 


where  x-.*  ,,,*  x are  the  quantities  (ranges  or  deflections)  in  the 
i n 

group  in  order  of  firing*  of  which  there  are  about  10*  and  i us  the 
arithmetic  mean  of  these  quantities.  Then  the  quantity 
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n~l 

Jl  * 1=1 

2(n  - 1) 

is  computed,  and  0.6745  min  ( a*  /v  is  accepted  as  being  suffi- 
ciently close  to  the  probable  error  with  the  effect  of  the  trends  in 
the  wind  and  velocity  removed 0*  Sometimes  the  above  quantities  are 
multiplied  by  a Beta  function  to  correct  for  sample  size. 

Footnotes  A basis  for  this  choice  is  the  following?  if  the  firings 

are  assumed  to  be  independent,  and  the  single  shot  expectations  for  g 

x_,  are  a.,  and  all  the  shots  are  assumed  to  have  the  same  variance,  <x“, 
i 1 

then  it  is  not  difficult  to  show  that  the  expected  value  of  Lq  is 


n«l 


*1  I*  a 2 +.  11 


^Vai+1) 


I in  particular,  for  a linear  trend  with 


2 (n«l) 


firings  at  equal  time  intervals,  a ® ai  * b,  we  would  have  that  the 

expected  value  of  Aq  is  and  that  of  is 

Now  ordinarily  a is  small  compared  to  a.  Hence,  in  this  case,  A^  is  a 
better  estimate  of  the  variance  than  Aq.  The  expected  value  of  the 
normalized  square  of  the  second  difference, 

n~2  n=“2  . 

*2  ■ <5i£37  ^ (vxi*2-2jW2»  (,i%2-2*i.i)2- 


Hence,  if  the  trend  is  linear,  this  expectation  is  v „ Were  it  not  for 

2 

the  small  number  of  rounds  (about  10 ) and  the  fact  that  in  practice  sT 

2 

is  a small  fraction  of  a , /L  would  be  a better  recommendation  for  an 

2 * 
estimate  of  c 0 


In  fact  if  the  trend  is  a polynomial  of  degree  k then  an  estimate 
using  k*l  differences  removes  the  trend  entirely,, 
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U,  Reduction  Trajectories 


Immediately  before  the  reduction  trajectories  are  begun  a first 
estimate  of  the  ballistic  coefficient  C is  made  from  the  average  mass 

for  each  group  and  the  caliber  of  the  projectile,  modified  if  experience 
dictates,  by  a form  factor  as  in  equation  (i;). 

A function  of  the  speed  of  sound,  a,  appropriate  for  use  in  either 

2 

(3)  or  (6),  e„g„5  1/a  , is  computed  by  hand  from  the  previously  reduced 
temperatures  using  (7)°  Using  this  function  together  with  the  previously 
reduced  p(y)  and  the  first  estimate  C^,  a machine  (Eniac,  Edvac,  or 

Ordvac)  computes  a reduction  trajectory  using  the  Heun  method  to  inte- 
grate  (l)» 

The  next  step  (see  the  block  flow  chart,  fig,  1)  is  to  make  a second 
estimate  such  that  the  range  x^  obtainable  with  it  will  lie  on  the 

other  side  of  the  observed  range  for  the  group  from  the  range  just  com- 
puted with  the  first  estimate  C^„  Then  another  reduction  trajectory  is 

computed  to  verify  or  improve  this  guess# 

It  has  been  found  from  experience  that  a sufficiently  close  brac- 
keting of  the  range  has  been  obtained  whenever  - 

(8)  -c2|  < .005  cz# 

Thus  if  (8)  is  not  satisfied  a new-:  estimate  of  is  ..made  by'.linear..  inter- 
polation (with  respect  to  the  ranges)  between  and  C2„  With  this  new 

C another  reduction  trajectory  is  computed  and  if  (8)  is  now  satisfied 
with  C2  and  the  new  estimate  then  we  go  on  to  another  set  of  reductions : 

for  another  pair  of  initial  elevation  and  velocity  values o If  (8)  Is 
not  satisfied  the  procedure  of  new  estimations  of  C and  computations  of 
reduction  trajectories  is  repeated  until  (8)  is  satisfied,,  The  number 
of  repetitions  seldom  exceeds  three. 

When  the  computation  bf  a set  of  trajectories  for  all  given  combi- 
nations of  muzzle  velocities  and  elevations  has  been  accomplished,  a 
fitting  of  the  CBs  finally  used  in  the  computations  by  a low  order 
polynomial  in  0 and  Vq  is  done  as  follows  0 First  a family  of  graphs 

of  C vs,  ff  with  as  a parameter  and  a family  of  graphs  of  C vs „ Vq 

with  0 as  a parameter  are  plotted.  An  experienced  person  determines 
visually  whether  a reasonable  fit  by  least  squares  can  be  made  by  a 

constant  C or  by  a linear  or  by  a quadratic  function  of  0 or  V . It 

o 

usually  occurs  that  0 appears  to  a higher  power  than  VQ  in  these  fitsj 
occasionally  Vq  does  not  appear  at  allj  furthermore,  the  power  of  0 
usually  does  not  exceed  20 


The  fit  is  a weighted  least  squares  fit?  given  by  minimizing 


C~C 


2 


P,E,obs 


».  V 


where  C(0,  V ) is  the  C obtained  from  the  reduction  trajectories,  C, 

■3=^  is  the  change  in  range  per  unit  change  in  C obtained  from  the  last 

few  reduction  trajectories  in  each  group,  and  P,  Egg  6cq)  is  the  ob~ 

served  probable  error  in  range,  The  criteria  for  the  degree  of  the  fit 
is  that  no  lower  degree  fit  produces  an  within  one  observed  probable 

error  of  the  computed  in  the  appropriate  reduction  trajectories « If 

it  so  happens  that  a quadratic  fit  is  not  obtainable  (as  may  happen  when 
high  angle  values  are  needed  in  the  table)  over  the  ertire  range  of  0 
and  V,  the  independent  fits  are  made  over  different  portions  of  the 
range  of  0, 


In  addition  to  the  quantity  dx.  /dC  the  quantities  cbc and 


dx^/30  are  obtained  from  the  last  few  reduction  trajectories.  They 


are  used  in  a statistical  formula  for  the  probable  error  in  range  due 
to  special  effects.  The  details  of  this  formula  are  too  specific  to 
be  discussed  here  beyond  remarking  that  it  involves  these  derivatives 
and  sanple  probable  errors  in  muzzle  speed,  mass,  elevation,  and  form 
factors.  If  these  special  effects  contribute  substantially  to  the 
observed  probable  error,  a new  jump  correction  is  made  to  reduce  the 
contribution  to  these  effects  and  a new  fit  to  C is  made.  Another  use 
of  these  quantities  is  made  later  in  obtaining  the  constants  in  a fitted 
formula  for  the  observed  probable  errors  in  range  and  deflection. 


E,  Normal  Trajectories  and  Differential  Effects 


Once  the  proper  ballistic  coefficient  is  found,  trajectories  are 
computed  under  standard  conditions  for  each  velocity  zone  and  a set  of 
angles  of  elevation..  These  standard  conditions  are  that  there  be  no 
Wind  (w^  * ^ * O),  no  rotation  of  earth  effect  (\^  « Vg  “ ^3  “ 0)9  no 

variation  in  g,  and  a standard  atmosphere.  At  present  a standard  atmos- 
phere, is  taken  as  one  for  which 


(9) 


=hy 

rs"  * 


p = p.  e * , a = ae  ",  h = ^ A, 


-Ay 


21 

o 


where  A is  a constant.  As  a result  of  these  conditions,  not  only  do 
the  first  equations  of  (l)  simplify,  but  alas  the  third  becomes  un- 
necessary since  the  trajectories  would  lie  in  a plane,  there  being  no 
sidewise  forces  and,  hence,  no  deflection. 
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These  normal  trajectories  give  range  and  time.cf  flight  for  each 
muzzle  speed  and  elevation*,  For  each  muzzle  speed,  interpolations  are 
performed  on  the  set  of  ranges  vs  elevations  to  give  range  as  the  in- 
dependent  variable  in  100  yard  intervals  with  elevation,  time  of  flight 
and  components  of  terminal  velocity  as  the  dependent  variable o 

From  this  interpolation  is  obtained  immediately  the  differential 
effect,  change  in  elevation  for  100  yard  change  in  range « 

Further  differential  effects  are  obtained  by  computing  a few 
additional  normal  trajectories  with  changes  in  air  density,  temperature, 
range  wind  and  muzzle  speed®  These  give  additional  items  for  the  tables 
and  for  later  use  in  computing  probable  errdrs  to  be  entered  in  the 
tables  and  to  be  used  in  computing  fork® 

During  the  course  of  the  computations  of  each  normal  trajectory 
the  values  of  x,  y,  and  y as  functions  of  t are  printed  out®  This 
is  to  be  used  later  in  obtaining  table  entries  for  angle  of  site  cor- 
rections  to  the  gun  elevation® 

One  also  gets  from  the  normal  trajectories  the  values  of  terminal 
speed  and  position  of  summit® 

F«  Firing  Table  Elements 

The  firing  tables  are  arranged  in  sections  with  a fixed  muzzle 
velocity  in  each  section®  Within  such  a section  the  independent  vari- 
able is  the  range,  given  in  hundred  yard  intervals®  We  have  already 
mentioned  how  the  corresponding  elevations  and  variations  in  range 
with  muzzle  velocity,  temperature,  wind  and  density  are  obtained  from 
the  effects  trajectories®  All  these  quantities  are  also  tabulated® 

Among  the  most  important  table  elements  remaining  are  corrected  time 
of  flight,  drift,  probable  errors  in  range  and  deflection,  fork,  angle 
of  fall,  and  complementary  angle  of  site  for  each  mil  of  site. 

The  corrected  time  of  flight  is  obtained  by  adding  to  the  computed 
time  of  flight  a correction  «^to  The  formula  for  this  correction  is 
gotten  by  a least  squares  fitting  of  linear  form  in  0 and  V (or  some 

other  simple  function  of  these  which  increases  more  rapidly  when  0 exceeds 

65°)  to  the  differences  between  the  observed  times  of  flight  and  the 
times  of  flight  of  the  reduction  trajectories® 

The  drift  is  the  portion  nf  the  deflection  which  is  not  due  to 
wind  or  rotation  of  earth  and  is  attributed  to  the  spin  of  the  pro- 
jectile. The  values  of  drift  in  the  table  are  obtained  from  a form- 
ula which  is  a least  squares  fit  to  the  drifts  in  the  form 

D a d t n 

63 
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or  in  the  form 

D = k tan  0 

where  k and  n are  fitted  constants,  n usually  being  near  2, 

The  effect  on  deflection  due  to  a unit  cross  wind  is  reasonably 
approximated  by  the  formula 

W„D„  - ~ v (t  - ^2  ) 


a> 


GS 


where  W„  D,  is  wind  deflection  angle,. 

The  formulas  for  probable  errors  in  range  and  deflection  are  now 
obtained.  The  probable  error  in  deflection  is  given  by 


° cos  0 


idle  re 


Xa> 


P » w f tan 

0 

where  (P is  the  weighted  observed  PoE,^  for  each  0O  The 
probable  errors  in  range  are  usually  given  by  a least  squares  fitting 


(PJS„„  f - (P.E.,r  f * (P.B.  r 


CO 


m 


! '±5i) 

£ C it?1 


rfAX»)  ^ + 9 a^Xcon  n o|  (p  y \2  /^xoo\2  _ \2  /AXo3^2 

J_(rr>  -rr  (P-Ey  (PJV  * 


There  appears  to  be  some  physical  basis  for  the  first  formula  above, 
but  none  for  the  second & Recently  the  Firing  Table  Branch  of  the  Com- 
nutbns  Laboratory  has  been  trying  to  fit  a drift  formula  (lull*  or  Lj.l5 
in  / ) dependent  on  the  speed  along  the  trajectory,  the  spin,  some 
physical  constants  of  the  shell  and  some  aerodynamics  coefficients  de- 
pendent on  Mach  number. 
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are  obtained  from  the  effects 


trajectories*  and  some  of  the  probable  errors  on  the  right  hand  side 
are  either  observed  or  are  obtained  from  the  reduction  trajectories  and 
the  others  are  fitted  constants  . 

The  range  probable  error  besides  being  a table  element  is  also  used 
to  obtain  fork,,  the  change  in  elevation  (in  mils)  for  a change  of  four 
probable  errors  in  range;  fork  is  given  by 


F 


lie  (PJ5„) 

100  < 


where  c is  the  change  in  elevation  for  100  yards  change  in  range « 


The  angle  of  fall  is  obtained  from  the  inverse  tangent  of  the 
ratio  of  the  vertical  and  horizontal  components  of  terminal  velocity 
obtained  in  the  computation  of  the  normal  trajectories,, 

The  complementary  angle  of  site  entries  in  the  table  are  obtained 
either  by  hand  or  more  recently  on  the  ENIAC  from  two  interpolations  on 
the  printed  values  of  x9  and  y in  the  normal  trajectory  computations „ 

The  first  interpolation  is  to  get  x as  the  independent  variable  in  in- 
tervals of  100  yards „ The  second  interpolation  is  on  the  values  of  y 
for  each  x to  obtain  from  which  one  gets  immediately  the  angle  of 

site  corrections o 

These  table  elements*  those  mentioned  in  the  previous  section* 
and  a few  others  not  mentioned  here  are  then  subtabulated  and  rearranged 
for  final  tabulation  and  typing  before  the  tables  are  reproduced  and 
bound o 


IV  - BOMBING  TABLES 

Ao  General  Bemarks 

Bombing  tables  are  bomb  aiming  data  used  by  bombardiers*  bomb  and 
sight  manufacturers  and  others « Paralleling  the  case  of  firing  tables, 
these  tables  are  produced  by 

a)  making  of  a series  of  bomb  drops  on  the  range  and  measuring 
whatever  is  necessary  including  non-standard  conditions  at  the 
time  of  range  drop* 

b)  computing  reduction  trajectories  to  obtain  a drag  coefficient 
Kjj.  and  a ballistic  coefficient  characteristic  of  the  particular 

bomb* 

c)  computing  a series  of  normal  trajectories  with  the  chosen 
and  ballistic  coefficient*  and 


d)  computing  various  correction  formulas  and  probable  errors  for 
use  in  corrections  for.  non-standard  conditions  occurring  at  time 
of  drop, 

A block  flow  chart  describing  these  steps  in  further  detail  appears 
in  Figure  20 

We  now  describe  what  is  done  in  the  main  blocks  of  this  flow  chart 
to  obtain  finally  the  bombing  table,, 


B0  Collection  of  Bange  Data 

At  the  bombing  range  groups  of  about  10  bomb  drops  are  made  for 
each  of  about  2 air  speeds  and  for  each  of  about  $ or  6 altitudes  of 
release®  Host  drops  are  made  from  airplanes  in  horizontal  flight® 

The  speed  of  the  airplane  at  bomb  released  and  the  release  coordinates 
are  obtained  from  ballistic  cameras  or  from  cameras  obacuraa,  the  choice 
depending  on  the  bombing  range  at  which  the  drops  are  made®  These  are 
synchronized  with  a chronograph  from  which  is  obtained  the  release 
time®  The  bomb  is  tracked  in  flight  by  means  of  Askania  cameras  also 
synchronized  with  the  chronograph®  The  Askania  films  furnish  azimuth 
and  elevation  of  the  bomb  vs®  time;  supplementary  information  as  a check 
on  release  time  and  coordinates  are  also  obtained  from  the  Askaniaso 
The  time  of  impact  is  obtained  from  geophones  synchronized  with  the 
chronograph®  From  Mitchell  or  Bowen-Knapp  cameras,  also  synchronized 
with  the  chronograph,  data  for  striking  velocity  and  supplementary 
data  for  time  of  flight  are  available®  The  impact  point  is  obtained 
from  theodolite  data® 

While  the  drops  are  going  on  several  meterological  data~gathering 
runs  are  made®  The  data  obtained  are  the  wind  components  down  range 
and  across  range  for  equal  intervals  in  time,  at  altitudes  up  to  the 
release  altitudes,  and  the  density,  temperature,  and  humidity  at  un- 
equal intervals  in  altitude  (equal  intervals  in  pressure)®  The  data 
are  obtained  as  on  the  firing  ranges,  using  instruments  in  balloons 
whose  ascensions  and  drifts  are  tracked  by  theodolites® 


These  ra*T  data  must  then  be  processed  by  preliminary  computations , 
film  readings,  and  reductions  to  obtain  functional  and  numerical  param- 
eters appearing  in  the  reduction  trajectories,,  The  purpose  of  the  re- 
duction trajectory  computations  here  is  not  merely  to  obtain  a ballis- 
tic coefficient  as  in  the  case  of  firing  tables  but  also  to  obtain  a drag 
coefficient  Kp(M)„ 


C0  Preliminary  Reductions 

The  numerical  parameters  are  the  mass  of  the  bomb  and  initial 
conditions  required  for  solving  the  trajectories  (l),  namely,  release 
coordinates,  xQ,  yQ,  2.q,  release  speed,  vq,  and  angle  of  release 

(usually  close  to  horizontal)* 

The  functional  parameters  required  ere  the  drag  coefficient  K_(M), 

the  air  density  p(y),  the  speed  of  sound  a(y),  the  absolute  temperature 

T(y),  and  the  wind  components  v (y),  V (y)»  The  speed  of  sound,  a(y), 

X z 

is  obtained  fromT(y)  as  in  equation  (7)0  The  density  p(y)  is  obtained 
from  T(y)  and  the  humidity « 

The  meteorological  data  obtained  are  reduced  and  interpolated  to 
give  air  density,  temperature,  and  range  winds  and  crosswinds  aloft  at 
equal  intervals  (500  ft„)  in  altitude,, 

The  various  films  of  the  release,  trajectory,  and  impact  are  read 
(on  Askania  film  readers,  Telereaders,  or  Iconologs),  the  resultant 
data  being  prepared  for  triangulation  and  interpolation  on  IBM  machines 
to  produce  on  OROTAC  position,  velocity,  and  acceleration  of  the  bomb 
at  various  time  intervals  from  release  to  impact,,  The  IBM  CPC  merely 
adds  tracking  corrections  to  the  readings  and  also  takes  differences 
from  these  results  to  check  for  film  reading  errors* 

The  positions,  velocities,  and  accelerations  obtained  are  then 
used  with  the  metro  data  in  the  trajectory  equations  to  get  K^(M)  at 

various  points  along  each  trajectory  from  the  equation 
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where  d is  the  diameter  of  the  bomb  and  m is  its  mass®  The  rotation 
of  earth  and  variation  in  g are  neglected  in  (10)  because  (lO)  is  to 
be  used  to  get  first  estimates  of  Kp(M)  to  be  used  in  the  reduction 

trajectories^  and  the  effects  of  these  phenomena  are  secondary® 

The  values  of  CM)  obtained  from  (lO)  for  each  drop  are  then 

plotted  against  M,  being  plotted  on  the  sane  paper  for  all  the  drops® 
Then  a curve  is  drawn  so  that  it  passes  close  to  the  averages  of  these 
points®  The  values  obtained  from  this  curve  are  then  the  values  of 
the  first  approximation  to  K^(M)  to  be  used  in  the  reduction  trajec- 
tories „ 

Reduction  Trajectories 

With  a choice  of  form  factor,  i,  equal  to  1 and  the  first  approx- 
imation to  just  described*  several  reduction  trajectories  from 

each  group  (with  similar  initial  conditions)  are  computed.  These  tra- 
jectories use  the  air  structure  observed  at  the  time  nearest  the  time 
of  the  drop  to  obtain  the  wind  forces  on  the  bomb  and  the  speed  of 
sound  a(y)j  they  also  use  the  variable  g given  by  equation  (2)® 

Next,  for  each  drop,  one  computes  at  the  same  altitudes  along  the 
trajectories  the  difference  between  the  corresponding  computed  and  ob- 
served ranges  and  times®  If  all  these  '‘residuals"  do  not  have  constant 
signs  along  the  trajectory,  the  curve  is  adjusted  until  all  residuals 

of  each  type  do  have  constant  signs  (usually  opposite  for  range  and 
time)®  Adjustments  of  the  curve  for  lower  values  of  M are  made 

using  data  from  trajectories  of  low  altitude  drops,  and  for  higher 
altitude  drops®  Then  the  form  factor,  i,  is  adjusted  by  either  an 
increase  or  a decrease  of  10 %,  depending  on  whether  the  computed 
ranges  are  over  or  short  with  respect  to  the  observed  ranges®  By 
linear  interpolation  a new  value  of  i i3obtained  which  yields  a 
better  approximation  to  the  observed  ranges® 

The  computed  ranges  arid  times  of  impact  for  each  group  are  com- 
pared with  the  observed  ranges  and  times  of  impact  by  comparing  the 
mean  residual  of  the  terminal  values  with  the  probable  errors  of  these 
same  residuals®  The  objective  is  to  get  each  of  these  bracketed  within 
the  probable  error  for  each  group®  If  any  group  mean  is  greater  than 
one  probable  errpr  for  the  group,  then  either  the  K^.  or  the  i is  ad- 
justed, The  adjustment  is  usually  made  to  reduce  as  much  as  possible 
the  mean  of  the  group  whose  mean  had  been  the  greatest  in  absolute 
value®  This  is  accomplished  by  changing  i,  which  is  done  by  using  the 
values  of  obtained  from  the  earlier  reduction  trajectories® 

If  most  of  the  mean  residuals  in  range  to  impact  axe  of  one  sign  and 
if  the  mean  residuals  in  time  of  impact  are  of  the  other,  adjustments 
in  i are  raadeo  Otherwise,  the  K_.ss  are  adjusted  as  before®  When 
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enough  bombing  range  data  are  provided  several  repetitions  of  the  above- 
described  procedure  suffice  to  attain  the  objective  mentioned  above „ 

The  last  value  of  i used  is  multiplied  by  the  laat  curve  used 
and  this  resulting  curve  is  defined  to  be  the  for  the  bomb*  Hence™ 
forth  the  form  factor  of  the  bomb  is  considered  to  be  1, 


E„  Normal  Trajectories  and  Differential  Effects 

Once  the  proper  for  the  bomb  is  founds  trajectories  are  computed 

under  standard  conditions  for  a large  number  of  altitudes  and  speeds 
of  release,  The  standard  conditions  are,  as  with  the  firing  tables,  a 
standard  atmosphere  as  in  (9),  no  wind,  and  no  rotation  of  earth  effects. 
However  g is  still  permitted  to  vary  with  altitude,  just  as  in  the  re~ 
duction  trajectories. 

Again  there  is  no  deflection,  and  hence  no  2-equation, 

These  normal  trajectories  give  range  and  time  of  flight  for  each 
altitude  and  time  air  speed.  Most  usually  there  is  no  variation  in  the 
angle  of  drop,  although  requests  for  "climb  and  glide"  information  do 
appear  from  time  to  time#  For  use  in  computing  the  table  elements  the 
terminal  values 9 x^,  x^3  y^  are  required  as  well  as  the  quantities 


range  lags 

B * T0  fl  / Vyc 

tine  lags 

^ " ^CO  f Vyco 

Differential  effects  due  to  wind  and  height  of  target  above  sea 

level  are  obtainable  from  normal  trajectories.  Information  for  the 

differential  effects  due  to  height  of  target -are  obtained  during  com^ 

putation  of  the  trajectory  by  printing  the  values  for  a height  of  £000 

feet.  The  differential  effects  due  to  wind  are  obtained  from  a five™ 

point  differentiation  formula  on  x vs,  v to  get  3I/3V,  How  this  is 

0)  o 

used  to  get  the  wind  effects  will  be  explained  later  in  the  sectioi*  on 
table  elements , 
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F0  Bombing  Table  Elements 

A bombing  table  contains  tabulations  of 

disk  speed  setting  (DS), 
trail  setting  in  mils  and  feet, 
time  of  flight  setting  in  seconds, 
and  their  differential  corrections 

all  as  variables  dependent  on 

true  air  speed  in  knots  at  10  knot  intervals, 
release  altitude  above  target  in  feet  at,  say,  100  foot  in™ 
tervals  at  low  altitudes  (e0g»$  below  12,000  feet)  and  at, 
say,  1000  foot  intervals  at  high  altitudes, 
and  the  ratio,  q,  of  differential  ballistic  wind  to  wind  at 
bombing  altitude,  at  intervals  of  0*1  from  0o0  to  0,7, 

Cq  * 0,1  is  omitted),, 

The  differential  corrections  are  due  to  the  effects  of  range  wind  and 
height  of  target  above  sea  level.  Also  included  is  a differential 
correction,  due  to  cross-wind  effects,  called  the  aiming  point  off- 
set upwind  in  feet  for  a 10  degree  drift.  There  is  also  an  entry- 
incidental  as  far  as  bombing  is  concerned  called  dropping  angle  and 
its  differential  corrections,  given  in  degrees,  as  functions  of 
ground  speed  and  release  altitude. 

At  present,  the  ENIAC  obtains  the  bomb  sight  settings  by  the 
following  procedures 

(a)  Use  a $ point  interpolation  on  the  values  of  and  t ob- 
tained in  the  computation  of  the  normal  trajectories  to  get  them  at 
1000  foot  intervals  of  release  altitude, 

(b)  Use  a £ point  interpolation  on  the  values  of  x and  t ob- 

(£>  c>> 

tained  in  (a)  to  get  them  at  regular  intervals  of  release  speed, 

dx 

(c)  Compute  and stb tabulate  values  of  k = t - ^OT  var^ous 

regular  intervals  in  release  altitude  and  release  speeS  and  the  two 

dr 

heights  of  target,  y^  * o and  * £000  feet  (g^00--  is  obtained  using 

o 

a five  point  formula  approximating  the  derivative,) 

(d)  Smooth  (ll  point  hth  order  formula)  t and  k with  respect  to 
VQ  to  obtain  them  at  the  release  altitudes,  release  speeds,  and  heights 
of  target  in  (c), 

(e)  Smooth  (as  in  (d))  t and  k obtained  In  (d)  with  respect  to 
release  altitude. 
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(f)  Compute  and  put  on  cards 


1.  tJ  » t ' “ k 9 (secs.) 

co  Co  q 

2.  DS J * 5300/t^!,  (disk  speed) 

3. 

lit  )l  e 1000/y  (mils),  where  y « y = y , 

^3 oq  ^ ^ 

5e  DS«  « 2.  (whenever  ^'p-  230  mils)# 

F x^«23y* 


6*  A«Pe0o  « *3,7365  X1  (ft*)  (Aiming  Point#  Offset)# 

7.  » tan*1  x Jy*  (dropping  angle)# 

j 4 W "t  ^ 

8*  DAg#  DA^  a tan”1  JSL^JLJL.f  *x  = 168*894  ft/sec# 


5300  (vn  + v ) 

— — ~ — — ~ST“ ~~  (where  v = + 20,  40,  60,  80, 

x * *23  y + w t 1 x 

<o  J x (a 

100  knots  in  ft/fiec,  whenever  ]ilp  230  mils)* 

(g)  For  bombing  table  elements  below  12,000  feet  from  values  of 
time  lag  A and  range  lag  B obtained  earlier  interpolations  as  in  (a) 
are  performed  to  obtain  A and  B in  100  foot  intervals  in  altitude j 
these  are  reconverted  to  time  of  flight  and  range  from  the  range  lag 
and  time  lag  equations  given  above* 


The  IBM  section,  at  present,  then  nerforms  the  necessary  subtrac- 
tion and  retabulations  to  get  data  in  form  for  printing  on  the  electro- 
matic  typewriter. 


0*  Probable  Errors 

Although  they  do  not  appear  in  the  bombing  tables,  probable  errors 
in  range  and  in  time  are  desired  for  the  general  report  on  each  program 
of  range  bombing  and  of  subsequent  production  of  bombing  tables*  There 
are  two  types  of  probable  errors  computed. 


One  set  of  probable  errors  in  range  and  in  time  is  obtained  from 
the  computation  of  the  last  reduction  trajectories*  with  the  accepted 
K_,  for  the  bomb*  and  is  used  as  a measure  of  the  goodness  of  the  fit 

of  the  The  probable  errors  are  computed  for  each  altitude  and 

speed  group.  They  are  obtained  from  the  sum  of  the  squared  deviations 
of  the  residuals  (range  or  time)  from  the  mean  deviation. 


The  other  set  is  a set  of  probable  errors  in  range  and  in  time  for 
each  group  of  bomb  drops  reduced  to  standard  conditions  and*  hence,  is 
obtained  from  extra  computations  of  normal  trajectories.  These  probable 
errors  are  used  as  a measure  of  the  goodness  of  the  performance  of  the 
bomb.  They  are  computed  as  follows „ From  the  residuals  used  in  com- 
puting the  first  set  of  probable  errors,  and  from  the  quantity  ^ ^ 

obtained  in  the  reduction  trajectory  program,  two  sets  (one  for  range 
and  the  other  for  time)  of  values  of  i are  obtained.  These  values  of 
i are  chosen  so  that  if  the  trajectories  had  been  computed  with  them 
instead  of  with  i = 1 the  residuals  would  have  vanished.  Now  extra 
normal  trajectories  are  run  with  values  of  i bracketing  the  above  to 
obtain  ranges  and  times . Efcr  interpolation  (linear)  ranges  and  times 
corresponding  to  the  above-mentioned  values  of  i are  computed.  The 
probable  errors  of  the  residuals  of  these  ranges  and  times  for. each 
group  are  then  computed. 


V - CONCLUDING  REMARKS 
A.  Comparison  of  Computational  Procedures 

It  was  remarked  at  the  end  of  the  introduction  that  the  block  flow 
charts  (Figures  1 and  2)  of  the  firing  and  bombing  table  computational 
procedures  reveal  basic  likenesses  and  differences.  The  mass  of  detail 
implied  in  each  block  of  the  charts  described  in  Sections  III  and  IV 
somewhat  obscures  these  comparisons.  It  is  of  interest  to  review  them. 

The  likenesses  appear  in  the  structure.  The  main  blocks  of  the 
charts  are  the  same  i data  are  gathered,  preliminary  reductions  are 
made,  reduction  trajectories  are  computed  to  reduce  to  standard  condi- 
tions, normal  trajectories  are  computed  at  standard  conditiohs  together 
with  differential  effects,  probable  errors  ate  computed,  and  then  the 
tables  are  produced.  The  chief  causes  of  the  likenesses  are  that  the 
physical  model  and  the  numerical  procedure  employed  are  essentially  the 
same  for  both.  These  were  described  in  Section  XI.  Finally  both  pro- 
duce terminal  values  under  standard  conditions  toward  the  end  of  the 
computations  (lie.  from  the  normal  trajectory). 

The  differences  appear  within  the  main  blocks  and  have  two  main 
causes.  These  are  the  measurement  methods  on  the  one  hand,  and  the 
different  kinematic  ccnditions  caused  by  the  different  tactical  re- 
quireraents  on  the  other. 
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As  has  been  briefly  remarked  in  the  Introduction,  in  ground  to 
ground  range  firings,  the  measuring  equipment  being  used  does  not  track 
the  missile©  Thus  only  initial  and  terminal  information  is  available 0 
Hence,  the  firing  table  computation  up  through  reduction  trajectories 
is  forced  to  be  a group  and  statistical  reduction  controlled  by  termi- 
nal conditions  and  producing  a ballistic  coefficient © The  remainder 
of  the  firing  table  computations  are  directed  toward  obtaining  tabu- 
lations of  parameters  and  differential  corrections  useful  to  the  man 
in  the  field  who  has  to  compute  gun  settings  (i0e©  the  azimuth  and 
elevation  for  firing)©  The  tactics  of  artillery  has  been  such  that 
it  has  not  been  necessary  to  do  these  computations  a few  seconds  be- 
fore firings  they  could  have  been  done  as  much  as  an  hour  before  c 
Furthermore  the  tactical  problem  is  such  that  the  whole  trajectory 
is  involved,  but  an  accurate  control  of  range  is  much  more  essential 
than  an  accurate  control  of  time  of  flight©  (This  should  not  imply 
that  an  accurate  knowledge  of  the  time  of  flight  is  not  desire able© 

It  is  needed  in  fusing)* 

On  the  other  hand,  in  the  air  to  ground  problem  of  bombing  from 
aircraft  the  measuring  equipment  tracks  the  missile*  This  permits  a 
round  by  round  reduction  producing  a when  the  reduction  trajectories 

are  completed©  However,  in  this  problem  only  the  down  part  of  the  tra- 
jectory is  involved,  but  the  time  of  flight  accuracy  is  just  as  essen- 
tial as  the  range  accuracy©  Also  the  tactics  of  bonding  require  that 
the  computation  of  the  release  time  be  done  within  seconds  preceding 
release 5 the  computing  of  release  time  must,  therefore,  be  automatic* 
Hence,  the  computing  done  after  the  normal  trajectories  are  finished 
is  directed  toward  tabulations  of  appropriate  settings  of  a bomb 
sights  which  automatically  gets  the  release  time* 

These  differences  account  for  the  different  processing  of  data 
and  the  different  timing  and  arrangement  of  personnel  to  do  the  work 
of  producing  the  tables*  The  difference  due  to  method  of  measurement 
can  easily  disappear*  As  an  example,  for  rocket  and  missile  compu- 
tations the  Firing  Tables  Branch  does  reductions©  The  possibility 

of  Doppler  and  accelerometer  measurements  has  already  been  mentioned 
in  the  Introduction*  However,  the  differences  caused  by  tactical  re- 
quirements can  be  much  more  essential*  As  an  example,  in  ground  to 
air  firing  (e*g»  anti-aircraft),  only  the  ascending  portion  of  the 
trajectory  is  involved,  and  time  of  flight  must  be  accurately  computed, 
whereas  range  is  of  lesser  importance,  especially  when  proximity  fuses 
are  used*  As  a last  example,  the  air  to  air  firing  tables  are  some- 
what different©  The  time  of  flight  is  of  utmost  importance,  range 
being  a minor  consideration©  Due  to  the  relatively  high  missile  speeds 
and  relatively  short  distances  between  gun  and  target,  in  this  case 
only  the  flatter  portions  of  the  trajectories^ are  computed©  Thus,  a 
Siacci  ballistic  theory  (see,  e.  g*.  Bliss  [l])  is  used  instead  of  the 
theory  mentioned  in  Section  II©  In  the  computation  of  some  of  these 
tables  much  more  than  in  other  types  of  tables,  the  problem  of  the 
large  yaws  encountered  by  ilssiles  fired  at  large  angles  to  the  di- 
rection of  the  plane's  motion  present  substantial  difficulties© 

This  introduces  another  variation  in  the  physical  theory  used  (see,  e©g*s 
Sterne [?]) © Finally,  the  table  elements,  slant  range  and  drop  needed 
here  ar  j lot  required  in  other  types  of  tables* 
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B „ Time  Requirements 

It  is  now  of  interest  to  have  a comparative  idea  of  the  amount  of 
time  required  in  each  of  the  blocks  of  figures  1 and  2. 

We  consider  the  firing  tables  first* 

The  range  firings  normally  take  from  three  to  six  weeks.  The  ini- 
tial reductions  at  the  range  usually  take  two  weeks  or  so,  but  overlap 
the  firing  time  except  for  maybe  a week  over.  The  further  reductions, 
before  rmning  reduction  trajectories,  consume  roughly  four  man  weeks 
for  a program  of  thirty  or  so  groups.  Such  a program  would  call  for 
about  one  hundred  and  fifty  reduction  trajectories.  While  each  such 
trajectory  takes  some  fifteen  minutes  on  the  EM  AC  (depending,  of  course, 
on  the  time  of  flight)  and  some  two  minutes  on  ORDVAC  (when  print-outs 
for  every  second  of  trajectory  are  not  required),  the  total  time  for 
the  reduction  trajectories  usually  consumes  from  one  to  five  weeks, 
depending  on  the  priority.  One  to  two  weeks  are  then  needed  to  fit 
for  ballistic  coefficients.  Some  three  thousand  normal  and  effects 
trajectories  are  then  needed,  computed  at  the  rate  of  about  80  to  100 
per  hour  on  ORDVAC  and  about  30  per  hour  on  ENIAC.  It  is  usual  to  do 
these  computations  over  a period  of  about  six  to  eight  weeks*  During 
this  time  probable  errors  and  incidental  fittings  are  computed.  The 
production  of  table  elements  then  requires  some  four  weeks.  Finally 
eight  to  ten  weeks  are  needed  to  transcribe  the  table  elements. 

This  completes  the  work  done  by  the  Computing  Laboratory, 

The  bombing  table  running  times  are  as  follows. 

The  range  bombings  and  appropriate  data  gathering  occupy  anywhere 
from  six  weeks  to  a year  and  a half  depending  on  the  size  and  urgency 
of  the  program,  with  many  taking  about  half  a year.  A small  program 
will  involve  about  thirty  drops  where  large  ones  have  around  two  hundred. 
It  takes  about  one  half  a man  day  to  do  one  ± ^auction  of  a ballistic 
camera  plate,  and  another  half  man  day  for  metro  data.  Two  to  three 
months  are  then  needed  to  read  the  Askaria  films. and  convert  the  readings 
to  the  form  from  which  the  initial  estimate  of  is  made*  This  con- 
version, described  by  (l0).,  takes  ab  out  ten  to  fifteen  seconds  per  point 
on  ORDVAC.  There  are  about  six  to  twelve  of  these  points  per  drop. 

This  amounts  to  a total  time  of  about  four  hours  or  so  of  computing. 

These  computations  are  done  at  scattered'  times  while  the  Askania  readings 
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are  going  ori„ . The  reduction  trajectories#  of  which  there  are  usually 
around  1*00  to  600#  take  less  than  four  days#  about  fifteen  being  done 
per  hour  on  ORDYAC  and  about  eight  to  ten  per  hour  on  ENT  AC  with  prints 
outs  every  second  of  trajectory.  Around  600  normal  trajectories  are 
needed « They  are  computed  at  about  100  to  130  per  hour  on  ORDYAC  and 
about  30  or  so  an  hour  on  ENIAC.  The  ORDYAC  computations  are  usually 
spread  out  over  about  two  to  three  shifts  while  on  ENIAC  they  are 
spread  ever  about  three  shifts*  The  probable  errors  take  about  fefur 
man  days  of  hand  computing  to  do  about  l5o  of  them  and  they  can  be 
done  while  the  table  elements  are  being  transcribed.  The  table  elements 
are  computed  on  ENIAC  in  about  fifty  to  sixty  hours.  The  transcription 
of  the  bombing  tables  is  usually  accomplished  with  about  one  day  of 
rearrangements  on  IBM  followed  by  about  two  weeks  on  the  electromatic 
typewriter*  When  this  is  done  by  hand  it  takes  about  one  and  one-half 
months  or  more  of  typing.  This  completes  the  work  of  the  Computing 
Laboratory. 


C . Improvements 

We  conclude  now  with  a sketch  of  improvements  recently  intro- 
duced or  under  active  consideration  and  further  suggestions  for  possible 
improvement  in  the  computational  procedures.  Unless  statement  is  made 
to  the  contrary#  all  suggestions  below  were  either  initiated  by  the  authors 
or  by  coders  within  the  Analysis  and  Computation  Branch  of  the  Computing 
Laboratory. 

First  of  all,  T„  E.  Sterne  suggested  in  August  1952  that  the^  process 
of  estimating  a first  in  bomb  reductions  be  simplified  by  not  using 

every  single  frame  of  an  Askania  film  record  of  a drop.  By  choosing  an 
unequally  spaced  set  of  frames  to  smooth  for  position,  velocity,  and 
acceleration  at  the  mid-point,  one  can  not  only  obtain  sufficient 

information  for  the  drop,  but  one  can  also  see  to  it  that  the  results 
at  different  points  are  independent  of  one  another,  thereby  eliminating 
possible  fictitious  oscillations  in  the  reduced  data.  The  choice 
suggested  reduces  the  number  of  Askania  readings  to  approximately  one 
sixth  of  what  it  was. 

Secondly#  the  physical  model  described  in  Section  II  is  being  modi- 
fied by  Mp.  H.  Reed  of  the  Computing  Laboratory^  Firing  Table  Branch  to 
a particle  trajectory  theory  that  includes  a lift  force  as  well  as  a 
drag  force.  This  is  being  tried  with  the  hope  that  in  the  case  of  high 
angle  trajectories  not  only  will  the  already  achieved  matching  to  time 
of  flight  and  range  be  obtained,  but  a better  fit  to  actual  height  of 
summit  will  result. 
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Thirdly,  whereas  hitherto  metro  data  for  bombing  have  been  inter- 
polated to  equally  spaced  altitude  intervals  as  a separate  computing 
operation  before  the  reduction  trajectories  are  begun,  the  reduction 
trajectory  computation  is  now  being  remechanized  to  work  directly  from 
the  unreduced  metro  data* 

Fourthly.,  another  mechanization  in  progress  is  concerned  with  the 
transcription  of  table  elements®  Up  to  now  the  table  elements  printed 
out  by  the  tabulators  attached  to  the  ENT AC  and  the  CPC  have  been  re- 
copied in  the  desired  format  by  a typist s taking  about  eight  to  ten 
weeks®  At  present  this  work  is  being  mechanized  in  such  a manner  that 
the  table  elements  are  so  arranged  on  the  output  cards  that  the  proper 
format  can  be  printed  out  directly  on  the  electromatic  typewriter  in 
a few  weeks®  (AH  that  is  l€ft  is  the  additional  typing  of  the  headings, 
introduction,  and  a few  small  hand  computed  tabulations®)  At  present 
this  mechanization  is  being  developed  on  the  ENIAC  and  can  possibly  be 
done  on  the  CPC® 

Fifthly,  in  connection  with  the  metro  data,  we  have  noted  the 

somewhat  more  complicated  metro  reduction  required  by  firing  table 

computation  because  of  the  fact  that  several  metro  runs  may  occur 

during  a single  group  of  firings®  It  is  suggested  that  nevertheless 

this  reduction  can  be  mechanized  without  much  difficulty®  Also  there 

should  be  no  difficulty  involved  in  mechanizing  the  transformation  of 

metro  temperature  data  to  a/a  . „ 

s 

Sixthly,  it  has  been  suggested  that  two  input  card  readers  instead 
of  one  could  make  the  ENIAC  and  the  IBM-CPC  much  more  efficient  and 
rapid  in  performing  interpolation  in  two  variables®  If  this  extra 
equipment  were  installed,  the  computation  of  the  initial  estimate  of 
the  Kp  for  bombs  obtained  from  the  metro  data  in  equal  altitude  inter- 
vals and  the  Askania  data  in  equal  time  intervals  would  be  done  much 
faster  on  the  ENIAC  and  could  even  be  done  on  the  CPC* 

Seventhly,  W„  W®  Leutert  has  pointed  out  that  with  the  intro- 
duction of  the  static  magnetic  memory  on  the  ENIAC  it  is  now  possible 
to  introduce  further  mechanizations  of  programs,  already  used  on  EDVAC 
and  CRDVAC,  onto  that  machine®  Two  examples  are  the  variation  of<£ht 
(see  end  of  Section  H)  and  the  bracketing  (see  beginning  of  Section 
HI  D)  of  the  ballistic  coefficient® 

Eighthly,  W,  Leutert  has  further  pointed  out  that  the 
G functions  in  the  firing  table  computations  be  approximated  by  ana- 
lytic expressions  for  use  on  the  machine  s;<J  In  this  connection,  see  the 
last  suggestion  below* 

Finally,  there  is  the  possibility  of  completely  mechanizing  the  1 
whole  procedure  of  both,  running  the  reduction  trajectories  and  deter- 
mining of  the  final  to  be  used  in  the  tables®  This  would  seem  to 

be  feasible  when  large  memories  become  available  on  the  machines® 
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With  the  present  memories  of  the  ERL  machines,  a partial  mechanization 
appears  possible  if  memory  requirements  for  the  are  reduced ® This 

would  seem  to  require  putting  in  either  only  a few  points  or  an  analytic 
representation  (e.,.go,  polynomial  segments)  for  the  curve  and  possibly 

doing  the  same  for  the  metro  datac  We  might  add  that  various  suggestions 
and  attempts  to  represent  and  G by  analytic  expressions  have  been 

frequently  made  and  for  a long  time®  (See,  e „go5  a suggestion  by  Mbrrey 

fig  ). 

At  present  the  following  improvement  suggested  by  W0  W®  Leutert  is 
being  attempted  for  bombing  table  reductions ® Each  Kq  curve  is  to  be 

approximated  by  several  low  order  polynomial  segments  using  least 
squares  fitting.,  The  coefficients  of  these  polynomials  (instead  of  the 
much  larger  number  of  tabular  values  of  the  Kjj.  curve)  are  then  to  be 

put  into  the  memory  of  the  machine „ This  permits  a program  for  the 
determination  of  i (See  figure  2 „),  all  the  raw  metro  data,  and  a pro- 
gram for  interpolation  in  the  retro  data  to  enter  the  memory  at  the 
same  time.,  This  will  save  the  present  separate  metro  reduction  runs, 
the  present  repeated  reading  in  of  the  metro  data  during  the  determi- 
nation of  is  and  the , punching  of  the  corresponding  's ® It  is  esti- 

mated  roughly  that  such  a procedure  will  save  an  hour  or  so  of  metro 
reduction  runs  and  about  one-third  to  one-half  the  machine  time  now 
spent  on  the  reduction  trajectory  program®  This  is  achieved  at  a cost 
of  no  more  than  about  two  man-days  of  least  square  fittings  for  each 

most  of  which  is  taken  up  by  the  preparation  of  input  data®  A further 
saving  is  in  the  considerable  reduction  of  the  possibilities  of  errors 
in  punchingof  input  data,  which  can  lead  frequently  to  wasted  machine 
time®  This  saving  is  difficult  to  measure  but  would  appear  not  to  be 
negligible®  Finally,  such  a method  opens  up  the  possibility  of  experi- 
menting to  see  whether  a sufficiently  accurate  cannot  be  obtained 

from  the  data  of  relatively  few  bomb  drops®  The  resulting  saving  would 
be  considerable,  since  dropping  bombs  from  aircraft  is  expensive® 
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ADDITIONAL  NOTEi  H OCO 

THE  BOXER  WITH  SHAOOW  CONTAIN  SIMILAR  \ 

OPERATIONS  TO  THOSE  IN  BOXES  WITH  ^ 

SHADOW  IN  BOMBING  TABLE  FLOW  CHART. 
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